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The bulk and shear viscosities (rj and () have been studied for quark-gluon-plasma produced in 
relativistic heavy ion collisions within semi-classical transport theory in a recently proposed quasi- 
particle model of (2+l)-flavor lattice QCD equation of state. These transport parameters have been 
found to be highly sensitive to the interactions present in hot QCD. Contributions to the transport 
coefficients from both the gluonic sector and the matter sector have been investigated. The matter 
sector is found to be significantly dominating over the gluonic sector, in the case of both rj and £. 
The temperature dependences of the quantities, C/S, and (/rj indicate a sharply rising trend for the 
( closer to QCD transition temperature. Both rj and £ are shown to be equally significant for the 
temperatures that are accessible in the relativistic heavy ion collision experiments, and hence play 
crucial role while investigating the properties of the quark-gluon plasma. 
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I. INTRODUCTION 

The study of transport coefficients for hot QCD matter 
is an area of intense research since the discovery of a 
fluid like picture of quark-gluon-plasma (QGP) in the 
relativistic heavy ion collider (RHIC) at BNL The 
discovery of the QGP is attributed to the fact that at 
extreme energy-density and temperature, the ordinary 
nuclear matter goes through a transition to the QGP 
phase as predicted by the finite temperature Quantum- 
Chromodynamics (QCD) (this transition is shown to be 
a crossover [2[ at the vanishing baryon density). 

To describe a fluid, shear and bulk viscosities rj and £ 
respectively) are very important physical quantities that 
characterize dissipative processes in its hydrodynamic 
evolution. The former describes the entropy production 
due to the transformation of the shape of hydrodynamic 
system at a constant volume, and the latter describes the 
entropy production at the constant rate of change of the 
volume of the system (hot fireball at the RHIC). More- 
over, tj/S, and C/5 serve as the inputs while studying 
the hydrodynamic evolution of the fluid [3, Q • One can 
also couple hydrodynamics with the Boltzmann descrip- 
tions at the later stages after the collisions of heavy ions 
at the RHIC, by maintaining the continuity of the entire 
stress energy tensor and currents. That could be trans- 
lated in terms of the viscous modifications to the ther- 
mal distributions functions of particles. This leads to a 
smooth transition from the hydrodynamic regime where 
the mean free paths are short to a region where hydrody- 
namics is inapplicable and Boltzmann treatments seems 
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to be justified Q . Therefore, this sets a way to study the 
impact of transport coefficients of the QGP in various 
processes at the RHIC, and the ongoing heavy ion ex- 
periments at Large Hadron Collider (LHC), CERN (e.g. 
dilepton production, quarkonia physics etc.). Regard- 
ing viscous corrections to dilepton production rate at the 
RHIC, we refer the reader to . The determination of r\ 
and £, has to be done separately from a microscopic the- 
ory, either from a transport equation Q with an appro- 
priate Force, Collision and Source terms or equivalcntly 
from the field theoretic approach using Grecn-Kubo for- 
mulas [|| (long wavelength behavior of the correlations 
among various components of the stress-energy tensor). 

The QGP is strongly interacting at the RHIC 
as inferred from the flow measurements, and strong jet 
quenching observed there. This observation is found to 
be consistent with the lattice simulations of the hot QCD 
equation of state(EOS) which predict a strongly 

interacting behavior even at temperatures which are a 
few T c (the QCD transition temperature) . The flow mea- 
surements suggest, a very tiny value for the ratio of rj 
to the entropy density, S(ri/S) for the QGP, and the 
near perfect fluid picture [Tlj-[l4j (except near the QCD 
transition temperature where Q/S is equally significant 

Ty/5 dMl). 

Preliminary studies at the LHC [i~9l42l1 | reconfirm 
above mentioned observations regarding the QGP. In 
heavy- ion collisions at the LHC, in addition to the el- 
liptic flow obtained at the RHIC, there are other inter- 
esting flow patterns, viz., the dipolar, and the triangular 
flow, which are sensitive to the initial collision geome- 
try [22|. There have been recent interesting studies to 
understand them at LHC [l]! [23| . A more precise mea- 
surements of various flows and jet quenching at LHC is 
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awaited. On the other hand, C has achieved considerable 
attention in the context of the QGP after the interesting 
reports on its rising value close to the QCD transition 
temperature [l5l |. Subsequently, the possible impact of 
the large bulk viscosity of the QGP at the RHIC have 
been studied by several authors; Song and Heinz (24| 
have studied the interplay of shear and bulk viscosities 
in the context of collective flow in heavy ion collisions. 
Their study revealed that one can not simply ignore the 
bulk viscosity while modeling the QGP. In this context, 
there are other interesting studies in the literature [25l — 
l3l| . The role of bulk viscosity in freeze out phenomenon 
has been offered in [2l], [32j . Effects of bulk viscosity 
in hadronic phase, and in the hadron emission has been 
studied in [331 ] . Interestingly, in the recent investigations, 
these transport coefficients are found to be very sensitive 
to the interactions [l3|, [l4[ , and nature of the phase tran- 
sition in QCD [34[ . Another crucial aspect of C is its 
influence on the domain of applicability of hydrodynam- 
ics at the RHIC, viz. the phenomenon cavitation. This 
phenomenon has been addressed in detail in the context 
of diverging value of C near the QCD transition temper- 
ature in [35l [36|. Thus, the determination of ?] and C 
for the QGP has multi-facet dimensions, and impact on 
the variety of physical phenomena at the RHIC and the 
LHC. 

The determination of 77 and C has been performed 
adopting the view-point based on the inference drawn 
from the experimental results, and the lattice QCD (the 
best known non-pcrturbativc technique to address the 
QGP). Lattice QCD has indeed been very successful to 
study the QGP thermodynamics. However, the the com- 
putation of transport coefficients in lattice QCD is a very 
non-trivial exercise, due to several uncertainties and in- 
adequacy in their determination. Despite that there are a 
few first results c omp uted from lattice QCD for bulk and 
shear viscosities [37H4Cf which have observed the small 
value of ij/S, and large CAS at the RHIC. A very re- 
cent interesting analysis [41| suggests that it is possible 
to compare the direct lattice results with the experiments 
at the RHIC. From such a comparison the QCD transi- 
tion temperature came out to be around 175MeV. More 
refined lattice studies on rj and C are awaited in the near 
future. 

The work presented in this paper is an attempt to 
achieve, (i) temperature dependence of r/ and C (The 
gluonic as well as the matter sector contributions to 
these transport parameters have been obtained by comb- 
ing a transport equation with a recently proposed quasi- 
particle model [HlH of (2+l)-navor lattice QCD EOS. 
Noteworthy point is that the matter sector has largely 
been ignored in the literature in this context), (ii) to 
understand the small r]/S, and large (/S for the QGP 
for the temperatures closer to T c . More precisely, in- 
puts has been taken from the computations of rj and C 
in quasi-particle models 0,0, [13, IS EH) an d combine 
the understanding with a transport theory determina- 
tion of them in the presence of Chromo-Weibel instabil- 



ities 0, S3, EH • The present work is the extension of 
our recent work on r\ 0, 0] , and C 0] for the gluonic 
sector, to the full QCD. 

The paper is organized as follows. In Sec. II, we 
present the formalism to compute the 77 and £■ The quasi- 
particle model and transport equation have also been dis- 
cussed in brief in the same section. In Sec. Ill, we have 
presented the results on the temperature dependence of rj 
and C in (2 + l)-flavor lattice QCD, and relevant physics. 
In Sec. IV, we have presented conclusions and future 
prospects of the work. 



II. DETERMINATION OF TRANSPORT 
COEFFICIENTS 

There may be a variety of physical phenomena which 
can lead to viscous effects in the QGP (or in general any 
interacting system) [B[. Among them, particular focus 
in on the viscous effects which get contributions from the 
the classical chromo-fields. 

The idea adopted here is based on the mechanism ear- 
lier proposed 0, 0] to explain the small viscosity of 
a weakly coupled, but expanding QGP. The mechanism 
in the context of the QGP is solely based on the particle 
transport processes in the turbulent plasmas 0] that are 
characterized by strongly excited random field modes in 
the certain regimes of instability. They coherently scatter 
the charged particles and thus reduce the rate of momen- 
tum transport. This eventually lead to the suppression of 
the transport coefficients in plasmas. This phenomenon 
has been studied both in electro-magnetic (EM) plas- 
mas [Hoj . and in non-abeli an p lasmas (QCD plasma) by 
Asakawa, Bass and Muller [L2l . l48l . and further employed 
for the realistic QGP EOS in 0[II|]. 

The condition for the spontaneous formation of tur- 
bulent fields can be achieved in EM plasmas with an 
anisotropic momentum distribution [5 II ] of charged par- 
ticles, and in the QGP with anisotropic distribution of 
thermal partons [52] ], In the context of pure SU(3) gauge 
theory, this mechanism turn out to be successful to ex- 
plain small shear viscosity of the QGP and larger bulk 
viscosity for the temperatures accessible at the RHIC and 
the LHC 0,0]. Here, extension has been desired to the 
case of realistic EOS for the QGP by incorporating the 
effects from the matter sector (quark-antiquarks). 

It will be seen later that the analysis leads to an in- 
teresting observation regarding the relative contribution 
of the gluonic and matter sector to the transport param- 
eters. Before, we give a brief description of the quasi- 
particle understanding of (2 + l)-lattice QCD that fur- 
nishes an appropriate modeling of equilibrium state. 



A. The quasi-particle description of hot QCD 

Quasi-particle description of hot QCD medium effects, 
is not a new concept. There has been several attempts 
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so far to understand hot QCD medium effects in terms 
of non-interacting/weakly interacting quasi-partons, viz., 
the effective thermal mass models [EH, [|54|, the effective 
mass model with temperature dependent bag parameter 
to cure the problem of thermodynamic inconsistency 54] , 
the effective quasi-particles with gluon condensate 55[ , 
the Polyakov loop models [56| (Polyakov loop acts as ef- 
fective fu gaci ty), and the quasi-partons with effective fu- 
gacities |42h44|. The last one, that will be employed 
here, shown to be fundamentally distinct from all other 
mentioned models, and in the spirit of Landau's theory of 
Fermi liquids. Moreover, the model has been highly suc- 
cessful, in interpreting the lattice QCD thermodynamics, 
and bulk and transport properties of hot QCD matter 
and the QGP in relativistic heavy ion collisions. 

In our quasi-particle description for (2+l)-flavor lat- 
tice QCD[4J], we start with the ansatz that the Lattice 
QCD EOS can be interpreted in terms of non-interacting 
quasi-partons having effective fugacities, z g ,z q which en- 
codes all the interaction effects, where z g denotes the ef- 
fective gluon fugacity, and z q , denotes the effective quark- 
fugacity respectively. In this approach the hot QCD 
medium is divided in to two sectors, viz., the effective 
gluonic sector, and the matter sector (light quark sec- 
tor, and strange quark sector). The former refers to the 
contribution of gluonic action to the pressure which also 
involves contributions from internal fermion lines. On 
the other hand, latter involve interactions among quark, 
anti-quarks, as well as their interactions with gluons. The 
ansatz can be translated to the form of the equilibrium 



distribution functions, 
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and /* as follows, 



fg 

•I eq 



fq 

J eg 



'"1 



z g exp(-/3p) 
1 - z g exp(-/3p) 

z q exp(-ffp) 
I + z q exp(—(3p) 



z q exp(— (3\/p 2 + TO 2 ) 



I + Z q CXp( — (3\/p 2 + TO 2 



(1) 



where to denotes the mass of the strange quark, which 
we choose to be O.lGeV. (3 = T _1 denotes inverse of the 
temperature. Here, we are working in the units where 
Boltzmann constant, Kb = 1, c = I, and h/2n = 1. We 
use the notation v g = 2(N 2 — 1) for gluonic degrees of 
freedom , v q — 2 x 2 x N c X 2 for light quarks, v s = 
2 x 2 x N c x I for the strange quark for SU(N C ). Here, 
we are dealing with SU(3), so N c = 3. Since the model 
is valid in the deconfmed phase of QCD (beyond T c ), 
therefore, the mass contributions of the light quarks can 
be neglected as compared to the temperature. Therefore, 
in our model, we only consider the mass for the strange 
quarks. 

The effective fugacity is not merely a temperature de- 
pendent parameter which encodes the hot QCD medium 



effects. It is very interesting and physically significant. 
The physical significance reflects in the modified disper- 
sion relation both in the gluonic and matter sector. In 
this description, the effective fugacities modify the single 
quasi-parton energy as follows, 



LU 9 p =p + T 2 d T ln(z g ) 
w| = V + T 2 d T ln{z q ) 
LJi = vV + m 2 + T 2 d T ln(z q ). 



(2) 



These dispersion relations can be explicated as follows. 
The single quasi-parton energy not only depends upon its 
momentum but also gets contribution from the collective 
excitations of the quasi-partons. The second term is like 
the gap in the energy-spectrum due to the presence of 
quasi-particle excitations. This makes the model more 
in the spirit of the Landau's theory of Fermi -liquids. 
A detailed discussions of the interpretation and physical 
significance of z g , and z q , we refer the reader to our re- 
cent work [44j]. Henceforth, we shall use gluon sector in 
the place of effective gluonic sector for the sake of ease. 
We shall now proceed to the determination r\ and £ in 
presence of Chromo-Weibel instabilities. 



B. Chromo-Weibel instability and the anomalous 
transport 

The determination of 77 and £ have been done in a 
multi-fold way. Firstly, we need an appropriate mod- 
eling of distribution function for the equilibrium state. 
Secondly, we need to set up an appropriate transport 
equation to determine the form of the perturbation to 
the distribution function. These two steps eventually de- 
termine these transport coefficients. For the former step, 
we employ the quasi-particle model for the (2+1 )-flavor 
lattice QCD EOS discussed earlier. 

The 77 and C has two contributions same as the shear 
viscosity in J12j , (i) from the Vlasov term which captures 
the long range component of the interactions, and (ii) the 
collision term, which models the short range component 
of the interaction. Here, we shall only concentrate on the 
former case. The determination of shear and bulk viscosi- 
ties from an appropriate collision term will be a matter of 
future investigations. Importantly, the analysis adopted 
here is based on weak coupling limit in QCD, therefore, 
the results are shown beyond 1.2T C assuming the validity 
of weak coupling results for the QGP there. Note that 
the interplay for anomalous and collisional components 
of r\ has been discussed in [T2T - fl4T ]. and in the case of 
£ for the pure gauge theory, a discussion has been pre- 
sented regarding the interplay of the collisional [57H59f , 
and anomalous components in fl8l ]. It seems that at the 
conceptual level all the observation in [l8| regarding the 
interplay will remain valid here. Since, we do not have 
results for the matter sector therefore, we shall not of- 
fer a quantitative discussions on such an inter-play here. 
There have been computations of transport parameters in 
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the case of pure gauge theory based on the effective mass 
models within relaxation time approximation [60| . The 
approach adopted, and the physics set-up is entirely dis- 
tinct in the present case. It is to be noted that the gluonic 
component in all the quantities is denoted by subscript g, 
light-quark components are denoted by subscript q, and 
strange quark component is denoted by subscript s 



C. Determination of £ and n 

Let us first briefly outline the standard procedure 
of determining transport coefficients in transport the- 
ory The bulk and shear viscosities, £ and 77 of the 
QGP in terms of equilibrium parton distribution func- 
tions are obtained by comparing the microscopic defini- 
tion of the stress tensor with the macroscopic definition 
of the viscous stress tensor. The microscopic definition 
of the stress tensor is as, 



d 3 P 
{2ir) d E p 



(3) 



On the other hand, macroscopic expression for the vis- 
cous stress tensor is given by, 



T lk = PS ik + eutiik - 2ri(Vu)ik - C,5 ik V ■ u, (4) 

where (V«)it is the traceless, symmetrized velocity 
gradient, and V ■ u is the divergence of the fluid velocity 
field. E p accounts for the dispersion relation. To deter- 
mine C an one writes the parton distribution functions 

as 



Zg/q exp(-/3u ■ p + r)) =F 1 



(5) 



Assuming that j\ (p, r) is a small perturbation to the 
equilibrium distribution, we expand f(p,r) and keep the 
linear order term in /1; this leads to, 



f(p,f f ) = fo(p)+Sf(p,r) 

= /o(p)fl + /i(p,f)(l±/o(p) ), 



where / = {f g ,f q ,f s }, and similarly /1 = 
{/fj/ii/il- The plus sign in the bracket is for gluons, 
and minus sign is for fermions (q and s). Next, we shall 
consider these quantities cxplicitlyin the gluonic and the 
matter sectors. As discussed in [l^, [HI, C and r\ arc de- 
termined by taking the following form of the perturbation 

A, 



Lo g T 2 > 



fi(p,r) = ~— -T^PiPj Ai 9 (p)Vu) y + A 2g (p)Si 



AW) 



PiPj ( A ls (p)Vu)ij + A 2s {p)Sij .(7) 



Here the dimcnsionless functions 

A.i g ,iq,is(p), A.2g t 2q,2s(p) measure the deviation from 
the equilibrium configuration. Ai(p), A2(p), lead to 77 
and C respectively. Note that A lgy i q ^ s (p) is a isotropic 
function of the momentum in contrast to A2 g ,2q,2s(p), 
which is a anisotropic in momentum p. 

Since £ and 77 are Lorentz scalars; they may be eval- 
uated conveniently in the local rest frame of the fluid 
(in the local rest frame fo = f eq )- Considering the 
a boost invariant longitudinal flow, V • u = £ and, 
(Vu)ij = j^diag(—l 7 —l 7 2) in the local rest frame. In 
this case, the perturbations, fi(p) take the form, 



flip) = 



Ai g (p) 

UJgT 2 T 



PI 



El 

3 



A 2g(P) p 2 



OJgT 2 T 

A 2q(p) 2 

A 2s(p) 2 



(8) 



where r is the proper time(r = \Jt 2 — z 2 ). The shear 
viscosities are obtained in terms of entirely unknown 
functions, A lgMM (p) as, 

% = T ^/05jA.„(ri«l + /.) 



Vq 



I5T 2 



' d3p -Pl 2 A lq (p)f q (l-f g ) 



57T° U>„ 



The bulk viscosities are obtained in terms of the un- 
known functions, A2g t 2q,2s(p), 





V 9 

3T 2 J 


f d 3 p 
8^ 


r> 2 

P 1 2 
— \P 




Vq 

3T 2 J 


r d 3 P 

8tt 3 


P\ 2 
— {P 




v s 

3T 2 J 


f d 3 P 
8tt 3 


r> 2 

P 1 2 
— \P 

UJ 2 



Notice that while obtaining the expression for the bulk 
viscosity, we have exploited the Landau-Lifshitz Condi- 
tion for the stress energy tensor. The factor (— 3c 2 -E 2 ) 
in the right-hand side of Ea. (|I0p is coming only because 
of that. For details, we refer the reader to 45j. Here, 



is the speed of sound square extracted from the lat- 
tice data on (2 + I)-flavor lattice QCD. The determina- 
tion of Ain .in .i.sfa), and 7y gi9iS can easily be done follow- 
ing [HEf, and A 2g.2q,2s{f) and Cg.,q,s following [U]. 
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D. Determination of the perturbative, Ai and A2 

To obtain a analytic expression for the perturbations, 
A12, in our analysis, one need to first set up the trans- 
port equation in the presence of turbulent color fields. 
This has been done in [l2T - [l4j in the recent past. Here, 
we only quote the linearized transport equation, with 
Vlasov-Duprcc diffusive term, which arise after consid- 
ering the ensemble average over turbulent color fields, in 
the light-cone frame. The transport equation thus ob- 
tained reads, 

V^f eq ( P ) + V A /l/e,(p)(l ± fe q ( P )) = 0, (f f ) 

where f eq = {f g ,f q ,f s }, and = p (the quasi-particle 
model does not change the group velocity of the quasi- 
partons). Note that Eq. (fTTj) is written in the absence of 
collision term, and assuming the weak coupling approxi- 
mation. 

The mathematical structure of the Vlasov-Dupree op- 
erator is as follows, 



V. = 



9 2 C 2 



2{N 2 - 1) U - 



< E z 



B 2 > T„ 



(12) 



where C2 is the quadratic Casimir invariant for partons. 
For gluons, C 2 = N c , and for quarks C 2 = (N 2 - l)/2N c . 
Here, oj = {ui g , u> q , u> s }, denotes the quasi-partons dis- 
persions. The parameter, r m is the time scale associated 
with instability in the field, and the operator L 2 is, 



L 2 = -(pxd p f 



(13) 

Since L 2 contains angular momentum operator L p , there- 
fore it gives non-vanishing contribution while operating 
on a anisotropic function of p. It will always lead to 
the vanishing contribution while operating on a isotropic 
function of p. Following [141 ] . the expression for the 



Ai g (p) is obtained as, 



2(iV 2 - 1W 



3C g g 2 < E 2 + B 2 > t„ 



(14) 



On the other hand, expressions for Ai gi i s are obtained 



as. 



Ai,(p) 



Ai.(p) 



2(N 2 l)u 2 q T 
3C /5 2 <E 2 + B 2 > T „ 
2(N 2 - 1)oj 2 T 

B 2 > T„ 



(15) 



3C /5 2 < E 2 

Now, we write the transport equation containing only 
those terms which contribute to bulk viscosity £ as, 

„2 



P 

3w 2 



(V-u)f eq (l±f eq ) 



T 
9 2 Q 



Following, [18(, we can obtain the mathematical forms 
of the corresponding perturbations, A 2 . We shall write 
down the expressions in the gluonic sector, and matter 
sector separately to avoid the confusion. The expression 
for A2 9 (p) is obtained as, 

4(N 2 l)Tc 2 p 2 2 PT 

= N c g 2 <E 2 + B 2 >r m p 2 ^^^ ^ )ln( 

(17) 

On the other hand, the expressions for A2 9 ,2s are ob- 
tained as, 

4(A/ 2 - l)Tu 2 (P 2 2 PT 
= C s g 2 <E 2 + B 2 >r m p 2 ^ ~ C ^ )H V^ ] 
A d 4(7V 2 - l)Tw 2 a ( p 2 2 2 p T 

^ M = C f g 2 <E 2 + B 2 >r mP 2 ^ - C ° ^ )H 7^ ] - 

(18) 

Here, C f = 

Next, we relate the denominator of Eqs. (fT"5"|) . p7|) . and 
(|18D to the parton energy loss parameter q = q gi q q , via 
the relation, 

Now the gluonic, contributions to rj, and £ in terms of 
q can be rewritten as follows, 

T 6 64 (N 2 - 1) 



c 9 = 



— PolyLog[6,z g ], 

q 67T Z 

4(iV c 2 - 1 ) ' ' 
3Tn 2 q 



J J p T dp T dp z {^ - c 2 w 2 ) 2 x 



ln(^)x/ 3 (l + / 9 ). 
Po 



(20) 



On the other hand, quark- antiquark viscosities in the 
matter sector are obtained as, 



64N 2 v q 



-^^ ){ - PolyL 9[ ,,- Zq]] 



64iV 2 i/ s 



3ir 2 q(N 2 - 1) 



PolyLog[6, —z q ] 



-—PolyLog[5, -z q ] 



N c v qtS 
3CfTir 2 q 

ln(^) x /,,,(] /,.,;. 
Pa 



1 7 i P 1 1 \i 
p T dp T dp z (— - c s u) q s ) x 



(21) 



The PolyLog functions that appear in the expressions for 
Vg,q,s are defined in terms of the series representation as, 



00 z- 
— X K 



Ploylog[n,x] = ^ — 



(22) 



fe=i 



3(7V 2 - 



< E 2 + B 2 > r m L 2 fi(p,r)f eq (l ± f eq fy hei:e n is a positive integer, and the convergence of 

the series is ensured by the fact that x < 1. Moreover, 
(l6folyLog[n, 1] = C(n), and also PolyLog[n, —1] ~ — C( n )- 
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III. TEMPERATURE DEPENDENCE OF 77 AND 

c 

The determination of 77, and £ in the gluonic and mat- 
ter sector, is incomplete unless to fix the temperature 
dependence of q in both the sectors. The determination 
of q has been presented in the various phenomenologi- 
cal studies [6l| . either based on the eikonal approxima- 
tion, or the higher twist approximation, at a particu- 
lar value of the temperature. Here, we choose the q for 
gluons as 4.5 Gev 2 / fm, and 2.0Gev 2 / fm for quarks, at 
T = OAGev [62| (this temperature, we denote as To). 
Since, q appears in the denominator in the expressions 
for i] and £• Therefore, any set of values higher then 
those mentioned above will further decrease the values 
of 77 and C- At T = T , we can see that q g = 2.2hq q . 
At this juncture, we do not know these parameters at 
all temperatures, so we assume this relation holds for all 
temperatures. This assumptions is based on definition of 
q in the leading order in hot QCD [63j . where its same 
for both gluons and quarks except that of the quadratic 
Casimir factor. We shall utilize the relation q g = 2.25q q , 
while studying the temperature dependence of various 
quantities in the next subsections The exact temperature 
dependence of q, employing the quasi-particle description 
of hot QCD is not known to us at the moment. This will 
be a matter of future investigations. 



A. Relative contributions 

In the section, discussions are mainly on, (i) relative 
contributions of various components of 77 with their ideal 
counter parts, (ii) gluonic verses matter sector for 77, and 

C- 

Note that the shear and bulk viscosities in the full QCD 
can be obtained by summing of all the individual contri- 
butions of the quasi-partons as, 

V = Vg + Vq+ Vs 

( = ( g + ( q +( s . (23) 

The additivity of various components here is attributed 
to the fact that all of them belong to same process viz. 
the anomalous transport. Viscosity contributions from 
distinct process (e.g. anomalous and collisional) are in- 
verse additive due to the fact that various rates [HI, |l8| 
are additive. 

Let us define the relative quantities of our interest. 
Firstly, the we shall define the ratios of various compo- 
nents of r\ to that for the ideal system of quarks and 
gluons (denoted as r) Id ), which are defined as follows, 



Rgi = 



Ra 



n Id 

Iq.s 



R, 



(Vg+Vg + Vs) 

Id\ 
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FIG. 1. (Color online) 77 relative to the that obtained using the 
ideal EOS for QGP, in the gluonic sector, and the full QCD 
is plotted as a function of T/T c . The solid curve denotes the 
gluonic sector and dashed line denotes the full QCD. Both 
R g i and Ri approach to the ideal limit asymptotically. 



Similarly, to compare the relative contributions among 
various components of 77, we define the following ratios, 

r = Ui 

J^gq — 



K gs= — 



R 



s< 1 



1q 

- Hi 

Vs 

Hi 



(25) 



On the other hand, to compare the relative contribu- 
tions among the various components of £, following quan- 
tities have been defined, 



R 9q i 
R gs i 
R sq = 



k 
k 



(26) 



(24) 



The quantities defined in Eqs. (|24j|26p have been 
shown as a functions of T/T Ci in Figs. [TJU The ra- 
tios Rgi and Ri are shown as a function of temperature 
in Fig. [T] The parameter q assumed to be same in the 
interacting and ideal sector. We have considered temper- 
ature dependence beyond 1.2T C . Both R g i, and Ri show 
that interactions significantly modify the shear viscosity 
in the gluonic sector and the full QCD at lower tempera- 
tures. Both of them lies within the range {0.40, 0.97} for 
the temperature range, {T/T c = 1.2,6.0}. R q i and R S i 
are shown in Fig. [2]as a function of temperature. Both of 
them sits at the top of each other. This is not surprising 
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servations suggest that rj could be thought of as a good 
diagnostic tool to distinguish various equations of state 
at RHIC and LHC. 



FIG. 2. (Color online)?? relative to the that obtained using 
the ideal EOS for the QGP, in the matter sector. The R q i is 
rj relative to r\ Id in the light-quark sector, and similarly R S i is 
for the strange quark sector. Both the curves sits on the top 
of each other since the mass effects from the strange-quark 
sector do not play significant role here. 



o 
o 

CD 
> 

CD 




FIG. 3. (Color online) Shear viscosity in effective gluonic 
sector relative to matter sector. The solid lines denotes r} g 
relative to rj q , thin dashed lines (middle) represents rj g relative 
to rjs, and upper thick dashed line represents rj a relative to r) q , 
as a function of T/T c . 



since the mass effects coming from the strange quark sec- 
tor contributes negligibly in the temperature range con- 
sidered here. The light quark sector and strange quarks 
differ with each other by a factor of two coming from 
the degrees of freedom. While considering the ratio R S i 
it cancels from the numerator and denominator. From 
FigEl h evident that the hot QCD interactions signif- 
icantly modify the shear viscosity in the matter sector 
same as in the gluonic sector as compared to the ideal 
counter part. All of them approaches asymptotically to 
the ideal limit which is nothing but 'unity'. These ob- 



Next, we investigate the gluonic shear and bulk viscosi- 
ties relative to that of the matter sector. The relevant 
quantities in this context of r\ are R gq , R gs , and R sq given 
in Eg. ((25)) . These are shown as a function of tempera- 
ture in Fig. El On the other hand for C, R qg , R qs , and 
R sq are shown as a function of temperature in Fig. 2] 
It can be observed from Fig. [3J and Fig. 0] that the 
matter sector contributions significantly dominates over 
the gluonic contributions as far as the 77 and C are con- 
cerned. This could perhaps be understood by the follow- 
ing facts, viz., the higher transport rates in the gluonic 
sector as compared to quark sector as encoded in q, and 
the interactions entering through the effective fugacities 
z g and z q . Quantitatively, rj g is ~ 0.125?7 g , and 0.250?7 S 
at T = 1.20T C , and increase quite slowly as a function of 
T/T c reaching around 0.135 r\ q around 6T C (see Fig. [3]). 
The r\ s almost stays 0.5 r\ q for the considered range of 
temperature (contribution from the strange quark mass 
is almost negligible). From Fig. |4j it can be observed 
that R 9q and R gs have same qualitative behavior as a 
function of temperature. The quantitative difference is 
because of a factor ~ 2, since £ s ~ 0.5£ g . Again the 
mass effects in the strange quark-sector play almost neg- 
ligible role. The ratio R 9q initially increases and attains 
a peak around T/T c ~ 1.37 and then decreases sharply 
until T/T c = 1.6 and slightly increases beyond 1.6 and in- 
dicating towards the saturation at higher temperatures. 
Quantitatively, Cg ~ 0.27£ g around 1.2T C , and 0.l3( q at 
3.0T C . These observations are very crucial in deciding the 
temperature dependence of rj and and the ratios rj/S, 
(/S and (/rj. Most of the recent studies devoted to the 
77 and C draw inferences for the QGP which are purely 
based on the study of pure SU(3) sector of QCD only. 
The matter sector has largely been ignored. In the light 
of the above observations, it is not desirable to exclude 
the matter sector since the dominant contributions are 
from there. Finally, we can obtain the exact value of the 
ratio rj/S, and (,/S by employing the values of q quoted 
earlier (q = A.hGeV 2 / fm for gluons and 2.QGeV 2 / fm for 
quarks at T = AOOMeV). rj/S thus obtained as 0.570 and 
C/5 came out to be 0.057 at T = 400Me^y As discussed 
earlier, to obtain the exact temperature dependence of 
77, and C, one requires to fix the temperature dependence 
of q within the quasi-particle model employed here. This 
will be taken up separately in the near future. The quan- 
tity which can be determined unambiguity is the ratio 
C/?7 which is very crucial in deciding when the hot QCD 
becomes conformal. In other words, until what value of 
the temperature the effects coming from £ are important 
while studying the QGP. We shall now proceed to discuss 
these issues next. 
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B. The ratio (/rj 

The behavior of the ratio £/rj as a function of tempera- 
ture is shown in Fig. [5j and the temperature dependence 
of the ratios ^ = tj§^ and ^ = is shown in Fig. [5] 

Most importantly, from Fig. [5J clear indications are 
observed that the £ in the gluonic sector, and full QCD 
diverge as we approach closer to T c (the results are not 
shown around T c , since such a quasi-particle picture may 
not be valid there.). The quantity C,/rj shows sharp de- 
crease until one reaches up to 1.4T C in the gluonic sector 
and 1.6T C in the full QCD sector. Beyond that the de- 
crease becomes slow and the ratio slowly approaches to 
zero. Such a behavior of C/?7 as a function of temperature 
could mainly be described in the formal expressions in 
Eqs. (f2"TJ)) . and (|2"Tj) , and decided only only by the temper- 
ature dependence of c 2 but also by the energy-dispersion 
relations, w g ,q, s , and the temperature dependence of the 
effective fugacities, z g ^ q . It is evident that there is no 
way to obtain a (c 2 — |) 2 factor out from the expres- 
sion while performing the integration. However, such a 
scaling could be realized whenever p « T 2 dT(ln(z g:q )), 
and io gi q lS happen to be independent of z g and z q , and 
the thermal distribution of quai-partons show near ideal 
behavior. It may perhaps be realized at a very high tem- 
perature which are not relevant to study the QGP in the 
RHIC and the LHC. Therefore, C,/r) obtained here does 
not follow either a quadratic scaling or a linear scaling 
with the conformal measure (c 2 — ^). The same con- 
clusions were obtained in the case of pure gauge the- 
ory recently [l8|. Note that for the scalar field theories, 
(/i] = 15(c 2 — |) 2 (quadratic scaling) and it has 

been found to be true for a photon gas coupled with the 
matter [65j . The quadratic scaling is also valid in the 
case of pcrturbative QCD with a proportionality factor 
different from 15 (66|. Furthermore, in the case of near 
conformal theories with gravity dual, C/77 shows linear 
dependence on (c 2 — |) [671 ]. 

Finally, in Fig. [6j ^ and ^ are plotted as a function of 
T/T c . For entropy density, we utilize the quasi-particle 
results which are shown to be consistent with the predic- 
tions of lattice QCD, and in all the plots c 2 has been 
obtained from the quasi-particle model employing the 
method quoted in |68j |. Interestingly, these are of same 
order at T = 1.2T C . Below that temperature the latter 
dominates over the former and vice versa for T > 1.2T C . 
The former increases, in contrast to latter as a function 
of T/T c . There is a sharp increase shown by the latter 
until one reaches 1.4T C , and beyond that the decrease 
is slower and one is quite close to the conformal limit of 
QCD. The important inference that could be drawn from 
here is that while studying the QGP one need to incor- 
porate the effects of both shear and bulk viscosity until 
approximately 1.5T C . This confirm our viewpoint that 
both 77 and £ have significant impact on the properties of 
the QGP at the RHIC and the LHC. 
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FIG. 4. (Color online) Temperature dependence of relative 
quasi-parton bulk viscosities. The solid line shows the behav- 
ior of ( g relative to £ g , thin dashed line shows the behavior of 
£ 3 relative to £ s , and the thick dashed line shows the behavior 
of £ 3 relative to ( q , as a function of T/T c . 
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FIG. 5. (Color online) Temperature dependence of the ratio 
(,/r), in the effective gluonic sector, and the full QCD. The 
solid line represents the gluonic sector, and the dashed line 
represents the full QCD. 



IV. CONCLUSIONS AND FUTURE 
PROSPECTS 

In conclusion, shear and bulk viscosities of hot QCD 
are estimated by combining a semi-classical transport 
equation with a quasi-particle realization of (2+l)-flavor 
lattice QCD. The effective gluonic sector contributes an 
order of magnitude lower as compared to the matter sec- 
tor while determining the transport coefficients of hot 
QCD and the QGP. This could perhaps be understood 
in terms of transport cross-sections of gluons and quark- 
antiquarks. Since transport coefficients are inversely pro- 
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FIG. 6. (Color online) The temperature dependence of the 
quantities, r)/S and (/S in the full QCD. Here 5* = ^ , 
where S denotes the entropy density. 



and £ came out to be highly sensitive to the presence of 
interactions. This can be visualized from the modulation 
of r), as compared to its ideal counter part, and large and 
rising value of £ for the temperatures that are closer to 
T c (due to large interaction measure there). 

The investigations on the other contribution to the 
shear and bulk viscosities (collisional etc.), and their in- 
terplay with the corresponding anomalous transport coef- 
ficients will be a matter of future investigations. It will be 
interesting to include the effects of non-vanishing baryon 
density to the transport coefficients of the QGP. More- 
over, one could include the anomalous transport coeffi- 
cients in the Boltzmann-transport theory approach and 
study the impact on the response functions and quarko- 
nia physics along the lines of [H, [7(| , as well as dilepton 
production at the RHIC and the LHC. These ideas will 
be studied in the near future. 
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